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1  .  Introduction 


Since  the  introduction  of  the  holographic  matched  filter^']  many  different  types  of  filters  for 
optical  correlators  have  been  proposed. P-41  xhe  matched  filter  is  designed  under  the  assumption 
that  a  wide  sense  stationary  input  noise  is  additive  (e.g.,  spatially  overlapping)  with  the  target.  The 
limitations  of  the  matched  filter  in  terms  of  broad  correlation  peak  and  sensitivity  to  distortion  are 
known.  P, 3]  Different  criteria  have  been  proposed  to  characterize  the  filter  performance  in  terms  of 
light  efficiency,  sharpness  of  the  correlation  peak,  and  peak-to-sidelobe  ratiot^’*’!,  and  the 
importance  of  finding  trade-offs  between  the  different  criteria  by  changing  the  filter  characteristics 
has  been  established. 

The  matched  filter  function  is  derived  under  the  condition  that  the  additive  input  noise  is 
overlapping  with  the  target  and  is  at  least  wide  sense  stationary.!- '3  in  many  pattern  recognition 
applications,  however,  the  input  scene  noise  does  not  overlap  the  target.  The  idea,  in  plain  words, 
is  that  the  target  is  in  the  foreground  and  blocks  the  scene  noise. 

In  this  report,  we  focus  on  the  development  of  optimal  correlation  filters  for  pattern  recognition 
with  spatially  disjoint  target  and  scene  noise.  We  analyze  the  consequences  of  this  assumption 
which  describes  many  realistic  physical  situations.  In  particular,  we  show  that  for  this  class  of 
problems  the  matched  filter  expressions  and  the  optimum  receivers  derived  under  the  overlapping 
input  signal  and  scene  noise  assumption  may  not  perform  well  in  the  presence  of  spatially  disjoint 
input  signal  and  scene  noise.  The  matched-filter  functions  are  derived  by  maximizing  the 
conventional  correlation  peak  signal-to-noise  ratio  (SNR)  metric.  Under  the  spatially  disjoint  target 
and  scene  noise  condition,  the  conventional  definition  of  SNR  metric  needs  to  be  modified  and 
new  optimum  receivers  and  filters  need  to  be  determined.  Owing  to  the  spatially  disjoint  input 
scene  noise  and  the  tai'get,  the  procedure  for  the  maximization  of  the  correlation  peak  SNR  will  not 
result  in  a  useful  optimum  matched-filter  function. t^l 

In  this  report,  we  present  four  techniques  to  detect  a  signal  in  the  presence  of  spatially  disjoint 
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scene  noise.  The  first  approach  uses  multiple  hypothesis  testing  to  describe  an  optimum  receiver  to 
detect  a  signal  in  spatially  disjoint  scene  noise.  The  second  approach  is  based  on  minimum  mean 
squai'e  error  filtering;  that  is,  a  Wiener  filter  is  designed  to  detect  a  signal  in  spatially  disjoint  scene 
noise.  The  third  technique  is  based  on  maximizing  a  performance  metric  to  design  an  optimum 
filter  to  detect  signals  in  spatially  disjoint  scene  noise.  The  new  performance  metric  is  defined  as 
the  ratio  of  the  output  signal  at  the  target  location  to  the  output  signal  energy  for  the  spatially 
disjoint  target  and  scene  noise  case.  The  fourth  approach  is  based  on  modifying  the  conventional 
definition  of  the  correlation  peak  SNR  to  take  into  account  the  non-overlapping  nature  of  the  target 
and  the  scene  noise.  A  generalized  matched  filter  is  designed  to  detect  a  signal  in  spatially  disjoint 
scene  noise  by  maximizing  the  modified  SNR  . 

In  section  2,  the  limitations  of  the  conventional  SNR  as  a  measure  of  the  matched  filter 
performance  for  spatially  disjoint  target  and  scene  noise  are  discussed.  A  limited  number  of 
computer  simulation  tests  will  be  provided  to  show  that  the  conventional  matched  filter  obtained  by 
maximizing  the  classic  definition  of  SNR  will  not  perform  well  for  this  case. 

In  section  3,  the  design  of  an  optimum  receiver  for  pattern  recognition  problems  with  input 
scene  noise  that  is  spatially  disjoint  with  the  target  is  described.  The  processor  is  designed  based 
on  multiple  alternative  hypothesis  testing. 

In  section  4,  a  minimum  mean  square  error  filter  for  pattern  recognition  problems  with  input 
scene  noise  that  is  spatially  disjoint  with  the  target  is  described.  The  filter  is  designed  to  have  an 
output  which  is  a  delta  function  located  at  the  position  of  the  target.  The  filter  minimizes  the  mean 
square  of  the  difference  between  the  desired  delta  function  and  the  filter  output  in  response  to  a 
target  in  the  presence  of  spatially  disjoint  scene  noise.  We  show  that  the  filter  output  has  a  well 
defined  peak  and  small  sidelobes  in  the  presence  of  disjoint  signal  and  scene  noise. 

In  section  5,  an  optimum  filter  is  designed  to  maximize  the  ratio  of  the  correlation  peak  to  the 
average  output  energy  for  a  signal  which  is  spatially  disjoint  with  the  input  scene  noise.  This  filter 
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takes  into  account  the  non-overlapping  nature  of  the  target  and  the  scene  noise,  the  effect  of  the 
non-whiteness  of  input  scene  noise,  the  limited  spatial  size  of  the  input  window,  and  the  unknown 
illumination  of  the  target.  We  will  show  that  the  filter  produces  a  shaip  output  signal  at  the  target 
location  with  a  low  output  noise  sidelobes.  We  show  that  this  optimum  filter  is  very  similar  to  the 
minimum  mean  square  error  filter  for  input  scene  noise  that  is  spatially  disjoint  with  the  target. 

In  section  6,  the  classic  definition  of  the  correlation  peak  SNR  is  modified  to  obtain  a  new 
metric  for  spatially  disjoint  target  and  scene  noise.The  modified  SNR  is  defined  as  the  ratio  of  the 
expected  value  squared  of  the  correlation  peak  to  the  spatial  average  of  the  output  variance  and  is 
denoted  SNR.  By  maximizing  SNR,  we  have  designed  a  generalized  matched  filter  function  that 
can  detect  a  target  in  the  presence  of  spatially  disjoint  scene  noise.  It  will  be  shown  that  the 
conventional  matched  filter  function  is  a  special  case  of  this  generalized  matched  filter  under  the 
condition  that  the  input  noise  is  overlapping  with  target,  and  has  infinite  size  and  zero  mean. 

2.  Limitations  of  the  Classic  Definition  of  the  Correlation  SNR  in  Pattern 
Recognition  with  Spatially  Disjoint  Target  and  Scene  Noise^^®^ 

In  this  section,  we  show  that  the  matched  filter  expressions  derived  under  the  overlapping  input 
signal  and  scene  noise  assumption  may  not  perform  well  in  the  presence  of  spatially  disjoint  input 
signal  and  scene  noise.  The  commonly  used  optimum  matched-filter  functions  are  derived  by 
maximizing  the  conventional  SNR.  We  discuss  the  limitation  of  the  conventional  SNR  as  a 
performance  measure  in  matched  filter  based  optical  pattern  recognition  for  input  scene  noise  that  is 
spatially  disjoint  (or  non-overlapping)  with  the  target.  For  matched  filters  the  classic  definition  of 
the  SNR  is  the  ratio  of  the  expected  value  squared  of  the  correlation  peak  to  the  variance  of  the 
peak.''^^  Under  the  assumption  that  the  input  scene  noise  is  stationaiy  and  overlaps  the  target,  the 
correlation  peak  variance  is  a  good  statistical  measure  of  the  variance  at  any  other  sample  point  of 
the  correlation  signal.  We  show  that  when  the  input  noise  is  spatially  disjoint  with  the  target,  the 
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correlation  peak  variance  may  not  be  a  good  measure  of  the  noise  fluctuations  over  the  entire  output 
plane.  For  clarity,  one  dimensional  notation  will  be  used  in  the  following  analysis. 

We  examine  the  correlation  SNR  generated  by  a  matched  filter  function  for  an  input  signal  s(t) 
that  contains  a  target  (reference  signal)  r(t)  in  the  presence  of  spatially  disjoint  additive  scene  noise 
n(t).  The  input  signal  and  its  Fourier  transform  are: 

s(t)=r(t)+n(t);  (j) 

IS((D)I  exp[7:t)s(m)]=  IR(cd)|  exp[/^j^((o)]  +  iN((o)l  expfji))j;;j((o)]  (2) 

Here,  IS(co)l  IR((o)l  expfj<))p,((o)],  and  IN((o)lexp|j  (co)j  are  the  Fourier 

transforms  of  the  input  signal,  the  reference  signal,  and  the  input  scene  noise,  respectively.  We 
assume  that  the  noise  n(t)  is  bounded  in  a  finite  window  and  that  the  Fourier  transform  of  its 
realization  exists. 

We  consider  nonlinearly  transformed  matched  filters. f^l  The  first  order  component  of 

the  k*  law  nonlinearly  transformed  matched  filter  that  generates  the  correlation  signal  is:f^J 


Hki(co)  =|R(co)|‘'exp{-j0R((o)},  (3) 

Here,  the  carrier  fiequency  is  dropped  for  simplicity,  and  k  is  the  severity  of  the  nonlinearity  (a 
given  constant).  A  family  of  filters  is  generated  by  changing  k.  A  conventional  matched  filter[l] 
can  be  produced  by  setting  k=l.  A  phase-only  filter[4,9]  is  produced  by  using  a  binary  nonlinearity 
(k=0). 

The  light  leaving  the  filter  plane  that  generates  the  first  order  correlation  term  is  the  product  of 
the  input  signal  Fourier  transform  and  the  filter  function  in  Eq.(3): 
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c  1  (co)=|R(cof  exp { -j(j)R(ff)) )  IR((l))lexp[j<t)j^(co)]+|R(co)|''exp{ -j(t)R(m) } IN((o)l  exp[/])N((o)]  (4) 


The  first  term  in  Eq.(4)  produces  the  correlation  term  and  the  second  term  produces  the  output 
noise  term.  It  is  assumed  that  carrier  frequency  is  large  enough  to  separate  the  higher  order  terms 
from  the  first  order  correlation  term. 

We  denote  the  filter  output  by  y(t).  The  conventional  correlation  peak  SNR  is  defined  as  the 
ratio  of  the  expected  value  squared  of  the  correlation  peak  amplitude  to  the  variance  of  the 
correlation  peak  amplitude: 


SNR= 


|E{y(0)}P 

Var{y(0)} 


(5) 


Here,  y(0)  is  the  correlation  peak  amplitude; 


|R((of+’dco+ 


|R((jL))HN((o)|exp  { j  [(1)n(co)-(1)r(co)]  }  dco  (6) 


In  Eq.  (6),  the  first  term  is  the  signal  term  and  is  deterministic.  The  variation  in  the  correlation 
peak  owing  to  the  noise  fluctuations  is  represented  by  the  second  term.  Using  Parseval's  theorem, 
the  contribution  of  the  noise  term  to  the  correlation  peak,  n^,  is: 


where  rp(t)  is  the  inverse  Fourier  transform  of  lR(co)|'^exp{j(t)R(co) ) .  Here,  u  is  a  union,  S  j  is  the  area 
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in  the  input  plane  that  covers  the  target  r(t),  and  S2  is  the  remaining  area  in  the  input  plane  that  covers 
the  input  noise.  The  second  equality  follows  since  n(t)  is  zxro  over  S  p  The  SNR  is 


SNR=- 


2tc 


|R(0)P'^'d(x)+  E[n(t)]iy(t)dt 


'(0eS2 


/ 

E{[  n(t)rk(t)dt-  E[n(t)]r|,(t)dt  ;P ) 


Tt)eS2 


i( 


(t)eS2 


(8) 


For  the  conventional  filter  (k=l),  we  have  r|(t)=r(t)  which  is  spatially  disjoint  with  the  noise 
term.  In  this  case,  the  contribution  of  the  noise  term  to  the  correlation  peak  is  zero;  that  is,  nQ  is 
zero.  Thus  the  denominator  of  the  expression  lor  Ihe  SNR  goes  to  zero.  Since  the  noise  does  not 
affect  the  correlation  peak,  the  correlation  peak  variance  is  zero,  and  the  correlation  SNR  of  the 
conventional  filter  is  infinite.  This  is  the  problem  with  the  classical  definition  of  the  SNR  under  the 
spatially  disjoint  condition  for  measuring  the  correlation  performance:  a  large  input  noise  could 
generate  large  output  noise  amplitude  at  points  other  than  the  correlation  peak  position  such  that 
false  alarms  aix  generated  indicating  an  overall  poor  system  performance.  However,  for  these 
cases  the  conventional  SNR  could  remain  high,  presenting  a  misleading  measure  of  system 
performance. 

For  nonlinearly  transformed  matched  filters  (kz^  I )  including  the  phase-only  filter,  rp(t)  "leaks" 
into  the  input  scene  noise  owing  to  the  Fourier  plane  nonlinear  transformation.  (Leaks  can  also 
occur  for  non-conventional  matched  filters. )  The  spreading  of  rp(t)  "grabs"  the  noise  and  "pulls" 
it  into  the  correlation  peak.  Let  us  consider  the  case  of  k  varying  from  one  to  zero.  As  the  severity 
of  the  nonlinearity  increases  (k  decreases),  the  degree  of  the  leakage  increases.  The  most  severe 
leakage  is  for  k=0  (phase-only  filter).  In  the  presence  of  leakage  the  noise  contribution  in  Eq.(7)  is 
the  product  of  rp(t)  and  n(t)  integrated  only  over  the  area  that  rp(t)  spreads  into  h(t).  This  means 
integration  over  the  leakage  area  within  S2.  However,  most  of  the  energy  of  rp(t)  may  remain  in 
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the  S I  area  which  means  that  the  portion  of  the  signal  in  the  leakage  area  may  be  small  and  that  the 
contribution  of  the  noise  fluctuation  to  the  correlation  peak  may  be  small.  As  a  result,  the  variance 
of  the  correlation  peak  would  remain  small  and  the  SNR  would  be  large  (although  naturally  smaller 
than  the  conventional  case,  which  is  infinity). 

Computer  simulations  are  performed  to  investigate  the  effects  of  the  input  scene  noise  on  the 
correlation  peak  SNR.  In  the  computer  simulations,  white  Gaussian  noise  is  added  to  a 
deterministic  reference  signal  (tank  image).  For  the  overlapping  case,  the  conventional  matched 
filter  is  optimum  for  wide  sense  stationary  white  Gaussian  noise.  The  additive  input  noise  is 
considered  for  both  the  overlapping  and  non-overlapping  target,  as  shown  in  Fig.  1.  The  expected 
value  of  the  noise,  mj^,  is  0.45.  The  standard  deviation  of  the  noise,  a^,  is  0.2.  Both  the 
conventional  filter  (k=l)  and  the  phase  only  filter  (k=0)  are  considered[please  see  Eq.(4)].  The 
phase  only  filter  results  in  the  most  severe  leakage  of  the  output  signal  into  output  noise.  In  the 
computer  simulation,  the  target  (a  tank  image)  and  the  input  image  are  sampled  to  reduced  size 
with  15x26  and  50x69  pixels  respectively.  The  size  of  the  input  array  is  128x256  pixels.  The 
reference  signal  amplitude  in  the  input  plane  is  normalized  to  have  a  maximum  value  of  unity.  The 
SNR  is  measured  for  the  first  order  correlation  peak  according  to  Eq.(8).  The  tests  were  repeated 
100  times  to  obtain  the  expected  values  and  to  measure  each  output  SNR. 

Using  the  target  in  the  spatially  disjoint  scene  noise  [please  see  Fig.  1(a)],  we  show  an  example 
of  the  output  correlation  signal  intensity  for  the  conventional  filter  in  Fig.  2(a).  In  the  tests,  we 
measured  the  conventional  SNR,  and  the  failure  frequency  in  the  detection  of  the  target.  The  failure 
frequency  indicates  the  number  of  times  that  any  pixel  value  of  the  output  noise  intensity  becomes 
larger  than  the  correlation  peak  intensity.  It  is  assumed  that  target  is  detected  successfully  when 
the  correlation  peak  is  larger  than  every  pixel  value  of  the  output  noise. 

For  the  noise  parameters  used  here  the  tests  performed  for  the  conventional  correlator  (k=l) 
produced  very  lai'ge  numerical  values  for  the  conventional  SNR  (infinite  in  theoiy  and  larger  than 
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10^0  in  computer  simulation)  but  failed  to  detect  the  tank  image  many  times  by  producing  a 
conelation  peak  to  sidelobe  ratio  of  less  than  unity.  The  failure  frequency  for  the  conventional  filter 
was  about  78  out  of  100.  On  the  other  hand,  the  kth  law  nonlinear  filter  for  k=0  (that  is,  the  phase- 
only  filter)  produced  a  much  lower  conventional  correlation  peak  SNR  (about  320),  but  it  was 
more  successful  in  detecting  the  target  (the  failure  frequency  was  zero).  It  is  evident  that  under  the 
spatially  disjoint  condition  a  large  correlation  peak  SNR  may  not  be  a  good  indication  of  the  system 
performance. 

We  have  peiformed  the  same  correlation  tests  with  the  same  target  and  the  same  scene  noise 
that  is  overlapping  with  the  target  [please  see  Fig.  1(b)].  The  correlation  signal  intensity  for  the 
conventional  filter  is  shown  in  Fig.  2(b).  In  the  tests,  the  conventional  filter  performed  very  well 
with  no  failure  and  produced  a  SNR  of  about  2.2x10-1  In  this  case,  the  input  noise  satisfies  the 
conditions  that  the  filter  is  designed  under,  and  depending  on  the  noise  color,  the  conventional 
expressions  for  the  matched  filter  result  in  optimum  performance.  The  effect  of  spatially  disjoint 
input  noise  mean  on  the  correlation  performance  is  discussed  in  Appendix  A. 

In  summary,  we  have  investigated  the  limitation  of  the  classic  definition  of  the  correlation  peak 
SNR  as  a  measure  of  system  performance  for  input  scene  noise  that  is  spatially  disjoint  (non¬ 
overlapping)  with  the  target.  The  correlation  peak  SNR  is  considered  for  the  classic  (conventional) 
matched  filter  and  for  the  nonlinearly  transformed  filters,  including  the  phase-only  filter.  The 
output  of  the  filter  contains  a  signal  term  and  a  noise  term.  For  spatially  disjoint  noise,  it  is  shown 
that  the  correlation  peak  SNR  is  infinity  for  the  conventional  filter  (k=l).  For  nonlinearly 
transformed  filters,  such  as  the  phase-only  filter,  the  output  signal  term  "leaks"  into  the  output 
noise  term.  The  degree  of  the  leakage  depends  on  the  filter  type.  However,  even  for  a  severe 
leakage  (phase  only  filter),  the  peak  SNR  remains  large  since  most  of  the  signal  energy  is  non¬ 
overlapping  with  the  noise.  Computer  simulations  that  use  spatially  disjoint  target  and  scene  noise 
are  provided.  The  filter  can  produce  a  large  numerical  correlation  peak  SNR  and  can  fail  to  detect 
the  target  by  generating  output  noise  larger  than  the  correlation  peak. 
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3.  Design  of  an  Optimum  Receiver  for  Pattern  Recognition  vrith  Spatially  Disjoint 
Target  and  Scene  Noised’*] 

In  Section  2,  it  was  shown  that  the  conventional  filters  designed  under  the  spatially  overlapping 
target  and  input  scene  noise  assumption  may  not  perform  well  for  input  scene  noise  that  is  spatially 
disjoint  (that  is,  non-overlapping)  with  the  target  J  We  know  from  classic  results  in 
communication  theory  that  using  hypothesis  testing  optimum  receivers  designed  to  detect  a  known 
signal  in  the  presence  of  additive  stationary  Gaussian  noise  can  be  implemented  using  matched 
filters.  For  pattern  recognition,  a  theoretical  design  of  an  optimum  processor  under  the  spatially 
disjoint  target  and  scene  noise  condition  is  needed.  In  this  section,  we  describe  the  design  of  an 
optimum  receiver  for  pattern  recognition  problems  when  the  input  signal  contains  a  target  in  the 
presence  of  input  scene  noise  such  that  the  input  scene  noise  is  disjoint  (non-overlapping)  with  the 
target  in  the  input  scene.  We  show  that  for  such  problems,  the  conventional  correlator  solution  that 
results  from  the  overlapping  input  signal  and  scene  noi,se  assumption  will  not  be  the  optimum 
solution.  We  shall  compare  the  optimum  receiver  obtained  for  the  spatially  disjoint  target  and 
scene  noise  case  with  the  classic  results  to  describe  the  additional  operations  that  need  to  be 
performed  on  a  matched  filter  to  obtain  the  optimum  receiver.  Since  in  general  detectors  produce 
an  additive  noise  on  the  input  image,  we  will  also  assume  that  the  target  can  be  noisy.  The  detector 
noise  is  added  both  to  the  target  and  to  the  background  noise.  However,  we  will  make  the 
assumption  that  the  noise  on  the  target  is  statistically  independent  of  the  noise  on  the  other  parts  of 
the  image.  This  assumption  is  easily  fulfilled  if  one  considers  that  the  detector  noise  is  white  and 
statistically  independent. 

We  use  multiple  hypothesis  testing  with  unknown  parameters  to  design  a  pattern  recognition 
system  for  the  spatially  disjoint  target  and  scene  noise  problem.  For  each  hypothesis,  the  received 
signal  s(t)  contains  a  target  (reference  signal)  r(t-tj)  at  location  tj  in  the  presence  of  spatially  disjoint 
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additive  input  scene  noise  nj(t).  Let  w,.(t)  denote  a  unit  magnitude  window  function  that  defines  the 
support  of  the  reference  signal  such  that  w^(.)  is  unity  within  the  reference  signal,  and  is  zero 
elsewhere.  It  is  evident  that  r(t-tj)Wj.(t-lj)=:  r(t-tj).  We  denote  the  detector  noise  by  n^(t)  and  we  assume 
that  it  is  statistically  independent  of  the  spatially  non-overlapping  background  noise  ng(t).  We  have: 

nj(t)=[nB(t)  +  nd(t)][l  -  w,.(t-tj)].  (9) 

If  the  detector  noise  is  negligible,  we  can  let  the  variance  of  n|.|(t)  go  to  zero,  and  the  target  will 
become  noise  free.  Since  the  amplitude  of  the  signal  r(t)  can  be  unknown,  r(t)  is  multiplied  by  an 
unknown  parameter  "a"  which  could  repre.scni  the  illumination  variation.  Different  hypotheses  H(j,a) 
correspond  to  a  possible  position  tj  of  the  target  and  to  a  possible  amplitude  parameter  "a": 

HCj,a):  s(t)=[ar(t-tj)+n^(t)]Wr(t-tj)  +  nj(t)  (10) 

This  is  an  estimation/detection  problem.  The  unknown  parameter  "a"  needs  to  be  determined 
followed  by  signal  detection.  We  use  the  Bayes  test  for  multiple  alternative  hypotheses  testing  to 
decide  which  hypothesis  is  more  probable,  and  to  determine  both  tj  and  "a".  Assuming  all  the 
hypotheses  equally  likely  and  equal  costs  for  all  errors,  the  situation  reduces  to  a  maximum 
likelihood  test,l'2]  and  the  decision  rule  becomes:  choose  the  hypothesis  H(j,a)  for  which  the 
corresponding  probability  density  function  (likelihood  function)  is  maximum.  The  assumption  of 
equally  likely  hypotheses  is  reasonable  since  we  do  not  have  the  probability  of  each  hypothesis.  We 
consider  "m"  amplitude  samples  of  the  observation  s(t),  where  s(ti)  is  the  observation  sample,  and  the 
vector  s  represents  the  m  samples.  We  denote  Sj  the  m^  samples  taken  over  the  target,  that  is,  within 
the  window  w,.(tptj).  We  denote  Tj  the  remaining  samples  taken  over  the  scene  noise,  that  is,  within 
the  scene  window  [l-w,.(tj-tj)].  The  joint  probability  density  function  of  the  m  samples  of  s(t)  given 


that  the  hypothesis  H(j,a)  is  true  is  denoted  by  p[s/H(j,a)].  We  denote  the  joint  probability  density 
functions  of  the  samples  of  n;(t)  and  nj(t),  as  Pn-(*)  and  Pn^jj-)-  respectively.  Since  it  is  assumed  that 
nj(t)  and  n|j(t)  over  the  target  are  statistically  independent,  p[s/H(i,a)]=pj|j(  s j)  pjij(sj-arj).  The 
vector  rj  represents  the  samples  of  the  target  r(tj-tj).  To  proceed  further,  the  probability  density 
function  of  the  detector  noise  Pnj(-)  is  needed.  We  shall  assume  that  njj(t)  is  white  Gaussian  with 
variance 

m  m 

Pnd(”d)=n  Pnd(nd(ti))=[(l/27tc>d^)®-^"']exp{-X  [nd(ti)]^/2c>d2}. 
i=l  i=l 

Under  these  assumptions,  the  log-likelihood  function,  log{p[s/H(i,a)] },  is: 

m 

log{p[s/HCj,5)]}=log[pnj(Sj)]-(l/2ad2)X  Wr(ti-tj)[s(ti)-Sr(ti-tj)]2  +  0.5mlog(]/27iad2).  (1 1) 

Note  that  in  Eq.  (1 1),  the  summation  in  the  second  term  is  over  the  entire  scene,  but  only 
samples  are  added  due  to  the  window  w^Ctj-tj).  To  choose  a  hypothesis  H(i,a),  we  set  to  maximize  its 
corresponding  log{p[s/H(j,a)]}.  The  last  term  of  Eq.  (1 1)  does  not  affect  the  maximization  and  will  be 
ignored.  The  maximum  likelihood  estimate  of  "a"  can  be  determined  by  maximizing  the  likelihood 
function  over  "a"  given  that  tj  is  fixed.  To  do  this,  we  set  to  zero  the  gradient  of  log{p[s/H(j,a)] } 
with  respect  to  "a"  and  obtain: 


m  X  s(ti)r(ti-tj) 

^  [s(ti)-ar(ti-tj)]r(ti-tj)=0.  or  a  =  - ,  (12) 

‘=’  I  r2(q-t,) 

i  =  l 


where  a  is  the  maximum  likelihood  estimate  of  "a"  given  tj.  We  substitute  the  maximum  likelihood 
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estimate  of  "a"  in  log{p[s/H(j,a)] )  of  Eq.  (11),  and  then  try  to  find  the  best  estimate  of  the  location  of 
the  target  by  maximizing  log{p[s/H(i,a)] )  over  tj.  Thus,  we  need  to  find  the  maximum  of: 

m 

log(p[s/H(j,a)]}=log[pr,.(Tj)]-(l/2ad2)^  w,fti-tj)[s(ti)-ar(ti-tj)]2  ,  (13) 

i=I 

where  log{p[s/H(j,a)] }  is  the  log  likelihood  lunction  using  a.  In  Eq.(13),  we  have  ignored  the  last 
term  of  Eq.  (1 1)  for  simplicity.  Expanding  the  right  hand  side  of  Eq.(13),  and  substituting  for  the 
estimated  value  of  "a"  given  by  Eq.(12)  in  Eq.(13),  the  optimal  decision  is  to  select  the  maximum  of 
the  following: 

log{p[s/H(j,a)]}=log[pn.(Sj)]  -  (l/2a(j2)Cjq(ti)  (14) 

where 

ni  2 

m  [Z  S(ti)r(ti-tj)J 

CN(tj)  =  Z  s2(ti)w,.(ti-tj)  -  -  . 

Z 

i=I 

Equation  (14)  provides  the  optimal  decision  rule  to  obtain  tj  and  needs  to  be  evaluated  for  all  tj. 

It  is  interesting  to  note  that  for  a  bounded  Piij(s j)  [that  is  0<pj^  .(Sj)<oo],  if  the  noise  on  the  target 
is  sufficiently  small  (o,^  «1),  then  the  actual  Pnj(Sj)  becomes  irrelevant  to  the  detection  process.  In 
this  case,  the  maximization  of  the  likelihood  function  over  tj  reduces  to  finding  the  maximum  of 
-CN(tj)  in  Eq.  (14).  Thus,  the  optimal  decision  to  obtain  tj  is  to  find  the  minimum  of  the  following: 
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(15) 


Ill  2 

„  f-X  .s(ti)r(ti-tj)] 

CN(tj)  =  X  s2(ti)w,-(ti-tj)  -  - 

'=‘  X  '■'(vtj) 

i=l 

When  the  input  scene  noise  is  spatially  disjoint  with  the  target  and  the  target  is  noise  free,  the 
optimum  receiver  is  independent  ol  the  scene  noise  probability  density  function.  In  Eq.  (15),  the  fiist 
term  computes  the  energy  ot  the  input  scene  over  the  w'indow  function  for  diffeient  tj.  The  second 
term  is  the  cross  correlation  betw'een  the  reference  function  and  the  input  scene  normalized  by  the 
reference  signal  energy.  This  normalization  is  due  to  the  presence  oi  the  unknown  amplitude 
illumination  "a".  The  minimum  possible  value  of  Cpj(tj)  is  zero,  which  value  occurs  only  when  tj  is  the. 
target's  true  location. 

As  an  aside,  since  C]s[(tj)>0,  Eq.  (15)  may  be  rearranged  as: 


[X  s(b)f(ti-tj)] 


M;(r)  = - - 

^  m  m 

X  I'^(ti-tj)X  s2(ti)w,.(ti-tj) 

i=l  i=l 


(16) 


where  there  is  again  equality  only  if  tj  is  the  target's  true  location.  In  general,  maximization  of  -Cfv[(tj) 
and  Mj(r)  will  not  yield  the  same  results;  however,  for  the  noise  free  target  which  is  disjoint  with 
scene  noise,  the  two  are  equivalent.  Equation  (16)  is  the  cross-correlation  between  the  target  and  the 
input  scene  normalized  by  the  energy  of  the  input  scene  over  the  window  function  Wj.(tptj)  and  the 
target  energy. H^l  It  is  evident  that  the  normalized  correlation  used  with  heuristic  arguments  is  optimal 
when  the  background  noise  is  spatially  disjoint  with  the  target  and  when  there  is  no  noise  on  the  target. 

It  is  interesting  to  note  that  if  the  noiseless  target  in  the  input  scene  is  a  known  to  have  an  amplitude 
scaling  of  one  (that  is,  a=l),  using  Eq.(l  1),  the  test  is  to  find  the  maximum  of: 
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(17) 


1 1 1 


This  test  implies  that  if  the  target  is  noise  free  and  is  known  exactly,  the  solution  is  to  find  t;  where 


1 1 1  1 1 1 

X  s(ti)r(ti-tj)  -O.sX  W,.(ti-tj)s2(q)  is  maximum.  This  occurs  for  t;  where  the  window  is  located  on  the 

i=i  i=i 

target  in  the  scene.  It  is  evident  that  the  conventional  correlation  is  not  optimal  in  this  case,  and  a  bias 

in 

sum  of  -0.52_^  w,.(tj-tj)s2(tj)  should  be  added  to  the  standard  correlation  receiver  X  ■'’(h)r(ti-t;). 
i=I  i  =  I 

If  the  noise  on  the  target  is  not  small,  then  the  probability  density  function  of  the  scene  noise 

Pnj(Sj)  is  needed,  and  the  normalized  correlation  will  not  be  optimal.  Let's  assume  that  Pnj(^j)  is 

Gaussian,  white,  with  mean  m^^  and  variance  ct]-|2,  and  that  (m-m|.)  statistically  independent  samples  of 

the  scene  noise  are  considered.  In  this  case,  the  decision  rule  is  to  find  the  maximum  of: 


.  ■■  »ii| 

log{p[s/Ha,a)]}=-[l/2(an2+a^2)]  ^  { [l-w,(ti-tj)],s(ti)-mn  )2  -  ( l/2ad2)CN(tj).  (18) 


A  limited  number  of  computer  simulation  tests  are  presented  to  investigate  the  performance  of  the 
optimum  receiver  designed  for  the  spatially  disjoint  input  signal  and  scene  noise.  The  tests  are 
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performed  to  find  the  location  of  the  target  according  to  the  optimal  decision  rule  presented  in  Eq.  (15). 
The  maximum  of  [-C]\[(tj)]  when  equal  to  zero  indicates  the  position  of  tj.  In  the  simulations,  white 
Gaussian  scene  noise  is  used  with  an  expected  value  of  m,-,  =  0.5  and  a  standard  deviation  of  Gf|=0.3. 
The  reference  signal  is  a  tank  image  as  shown  in  Fig.  3.  The  size  of  the  target  (tank  image)  is  15x26 
pixels,  the  size  of  the  input  signal  is  50x60  pixels,  and  the  size  of  the  input  array  is  128x128  pixels. 
The  reference  signal  amplitude  in  the  input  plane  is  normalized  to  have  a  maximum  value  of  unity.  The 
reference  signal  energy  is  28.3  and  its  mean  is  0. 14.  An  example  of  the  noise  free  target  with  non¬ 
overlapping  scene  noise,  is  shown  in  Fig.  3(a).  The  target  has  an  illumination  ol  0.5  which  is 
unknown  to  the  receiver.  For  this  case,  we  show  an  example  of  [-Ci\[(tj)]  in  Fig.  3(b).  The  maximum 
value  is  equal  to  zero  and  indicates  the  location  of  the  target.  The  test  was  repeated  many  times  and  the 
target  detection  was  always  successful  (as  it  should  be).  It  is  assumed  that  target  is  detected 
successfully  when  the  position  of  the  maximum  of  [-Cjqftj)]  when  equal  to  zero  coincides  with  the 
actual  location  of  the  target.  In  addition,  in  Fig.  3(c),  we  have  shown  the  cross-correlation  between  the 

m 

target  and  the  input  scene,  that  is,  ^  s(t)r(t-tj)dt,  for  the  disjoint  signal  and  scene  noise  of  Fig.  3(a). 

i=i 

Note  that  this  is  the  optimum  receiver  solution  under  overlapping  signal  and  scene  noise  assumption 
for  additive  white  zero  mean  Gaussian  noise.  It  is  evident  from  this  figure  that  the  cross-correlation  test 
by  itself  may  not  be  a  good  test  to  determine  the  location  and  presence  of  the  target  for  the  disjoint 
signal  and  scene  noise  problem. 

In  summary,  we  have  used  hypothesis  testing  to  design  an  optimum  receiver  for  pattern 
recognition  problems  with  input  scene  noise  that  is  spatially  disjoint  (or  non-overlapping)  with  a  noisy 
target.  We  have  used  white  Gaussian  statistics  for  the  detector  noise  to  illustrate  the  optimum  receiver. 
The  cases  of  a  noise  free  target  which  is  spatially  disjoint  with  the  scene  noise  and  of  a  target  which 
has  an  unknown  illumination  scale  are  discussed.  It  is  shown  that  for  a  noise  free  target,  the  optimum 
receiver  is  similar  to  a  correlator  normalized  by  the  input  scene  energy  within  the  target  window.  In 
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addition,  when  the  detector  noise  is  sufficiently  small  (noise  free  target),  and  given  that  the  scene  noise 
probability  density  function  is  bounded,  then  the  actual  scene  noise  statistics  becomes  irrelevant  to  the 
detection  piocess.  These  lesults  clearly  indicate  that  the  spatially  non-overlapping  target  and  scene 
noise  case  lesults  in  optimum  solutions  that  ai'e  very  different  from  the  conventional  solutions  obtained 
by  the  additive  spatially  overlapping  signal  and  noise  case. 

4.  Minimum  Mean  Square  Error  Filter  for  Pattern  Recognition  with  Spatially 
Disjoint  Target  and  Scene  Noiset'-^  '^^l 

In  this  section,  we  describe  a  minimum  mean  square  error  filter  for  pattern  recognition 
problems  with  input  scene  noise  that  is  spatially  disjoint  (or  non-overlapping)  with  the  target.  To 
piovide  a  sharp  output  lesponse  and  a  quiet  background,  the  filter  is  designed  to  have  an  output 
which  is  a  delta  function  located  at  the  position  of  the  target  in  the  input  scene.  The  filter  minimizes 
the  mean  square  of  the  difference  between  the  desired  delta  function  and  the  filter  output  in 
response  to  a  noisy  input.  Under  the  disjoint  signal  and  scene  noise  condition,  the  minimum  mean 
square  error  filter  solution  will  be  dependent  on  the  window  function  that  defines  the  target  and  it 
will  be  different  from  the  conventional  matched  filter  function  and  the  Wiener  filter.  Computer 
simulation  is  provided  to  show  the  performance  of  the  minimum  mean  square  error  filter  in  the 
presence  of  the  spatially  disjoint  target  and  scene  noise.  We  show  that  the  filter  output  has  a  well 
defined  peak  and  small  sidelobes  in  the  presence  of  disjoint  signal  and  scene  noise. 

In  the  analysis,  samples  of  the  signals  are  represented  in  discrete  form  and  are  used  in  the 
analysis.  Samples  of  the  input  signal  sequence  is  s(t|)  where  i=...,-l,  0,  1,  ...  is  an  integer.  It  is 
assumed  that  the  input  scene  has  a  infinite  size,  that  is,  the  input  signal  sequence  s(i)  is  infinitely 
long.  The  reference  signal  sequence  r(tj)  in  the  scene  has  a  finite  size  and  has  a  window  wr(ti).  The 
window  sequence  Wr(tj)  of  the  target  r(ti)  is  defined  to  be  one  wherever  r(tj)  exists  and  to  be  zero 
otherwise.  The  non-overlapping  scene  noise  is  defined  as  n(tj)[l-wr(ti)],  and  the  noise  n(ti)  is 
assumed  to  be  samples  of  a  wide  sense  stationary  random  process.  The  target  is  assumed  to  be 
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located  at  the  origin  in  the  input  scene.  The  input  scene  s(tj)  is: 


s(ti)  =  r(ti)+  n(ti)  [1  -  Wr(ti)]  =  r(ti)+  n(ti)  w(tO 


where  w(ti)=[l  -  w,-(tj)] 


Let  us  denote  the  impulse  response  of  the  minimum  mean  square  error  filter  by  the  sequence 
h(ti),  the  actual  filter  output  response  by  y(tj)  =  s(tj)*h(ti),  and  the  ideal  output  response  by 
y(tj)  =  5(i).  Here,  is  the  convolution  operation  and  5(i)  is  the  Kronecker  delta  sequence.  We 
use  the  minimum  mean  square  error  criterion  to  optimize  h(tj).  The  mean  square  error  (MSE)  is  the 
error  squared  between  the  ideal  output  and  the  actual  filter  output. 

MSE  =  X  MSE(i)  =  X  E{ly(ti)  -  y(ti)p} 

;  i 

where  MSE(i)=|y(ti)  -  y(ti)p.  In  this  section,  we  consider  that  the  signal  is  periodic 
(cyclostationaryf^^)  and  we  sum  over  one  period.  This  allows  to  avoid  divergence  problems  in  the 
MSE  dertermination.  Here  I  .  1  denotes  absolute  value,  and  E{.}  the  expected  value.  From  Eq. 
(20),  it  can  be  shown  that  the  mean  square  error  MSE  is: 

MSE  =  1+X  X  h(ti)b(ti,tj)h(tj)-2X  h(ti)E{s(ti)}  (21) 


where 


b(ti,tj)  =  X  E(s(ti+tp)s(tj+tp)) 

p 

=  Rrr(ti-tj)  +  mnRwr(li-Ij)  +  mnRrw(h-Ij)  +  Rww(ti-tj)Rnn(li'Ij)  =b(ti-tj) 
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Here,  R‘^(.)  is  the  correlation  function.  The  subscripts  rr,  wr  and  ww  denote  autocorrelation  of  r(.), 
cross-correlations  between  w(.)  and  r(.),  and  autocorrelation  of  w(.)  respectively. 

Using  Eq.  (21),  the  mean  square  error  MSE  can  be  written  as: 

MSE=  1+^  H(v)B(v)H''(v)-22^  H(v)E{S(v)},  (23) 

V  V 


where, 

B(v)=  |R(v)p+  mnW(v)R‘'(v)  +  mnW*(v)R(v)  -i-  |W(v)p*Sn(v)  (24) 

where  W(v)=[5(v)-mnWr(v)]  is  the  discrete  Fourier  transform  of  w(tj).  Here,  H(v),  B(v),  S(v), 
R(v)  and  Wr(v)  are  discrete  Fourier  transforms  of  h(tj),  b(tj),  s(ti),  r(tj)  and  Wi-(ti),  respectively. 
Sn(v)  is  the  power  spectral  density  of  the  input  noise  sequence  n(tj). 

The  MMSE  filter  is  obtained  by  minimizing  the  MSE  [plea.se  see  Eq.  (23)]: 

H‘(v)=S.^W> 

B(v) 


_ _  R(v)  +  mnW(v) 

|R(v)p-t-  mnW(v)R*(v)  +  mnW*(v)R(v)  4-  |W(v)p*Sn(v) 

The  input  noise  sequence  n(tj)  can  be  considered  as  a  sum  of  a  zero  mean  noise  sequence 
no(ti)  and  its  mean  value  m^,  that  is,  n(tj)=mj,-t-nQ(tj).  The  non-overlapping  input  scene  noise 
n(ti)w(ti)  can  thus  be  decomposed  into  two  parts,  m„w(tj)  and  n^)(tj)w(ti).  The  first  component 

deterministic  when  the  mean  of  the  noise  is  known  and  the  uncertainty  comes  from  the 
second  component  nQ(ti)w(ti).  The  input  data  can  be  written  as: 
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s(tj)  =  {r(tj)+mnW(ti)  )  +  {nQ(ti)w(ti)} 


(26) 


The  first  term,  {r(ti)+mnW(tj)},  can  be  considered  as  the  signal  to  be  detected  while  the  second 
term  {nQ(ti)w(tj)}  is  the  noise  term  to  be  rejected.  It  is  evident  from  Eq.(25)  that  the  numerator  of 
the  filter  transfer  function,  R(v)  +  mnW(v),  is  matched  to  detect  {r(tj)+mnW(ti)).  When  the  scene 
noise  is  spatially  disjoint  with  the  signal,  the  window  of  the  reference  signal  becomes  as  important 
as  the  reference  signal  itself.  This  is  extremely  useful  when  the  mean  of  the  noise  background  is 
strong.  In  this  case,  the  disjoint  window  {w(ti)=[l-Wj.(tj)])  may  contribute  significantly  to  the 
output  signal. 

In  correlation  tests  using  spatially  disjoint  target  and  scene  noise,  we  have  observed  that  when 
the  scene  noise  has  large  mean,  most  of  the  widely  used  filters  designed  lor  overlapping  signal  and 
scene  noise  may  fail  to  detect  the  target.  These  filters  do  not  utilize  the  di.sjoint  window 
information  of  the  reference  signal  and  scene  noise.  On  the  other  hand,  our  tests  indicate  that  the 
minimum  mean  square  error  filter  presented  here  works  very  well  for  spatially  disjoint  target  and 
scene  noise. 

The  image  in  Figure  4(a)  is  used  as  the  input  data  in  the  computer  simulation.  The  reference 
signal  is  a  tank  image  of  15x26  pixels  and  normalized  to  a  maximum  of  unity  as  shown  in  Fig. 

4(a).  The  target  (tank  image)  is  placed  at  pixel  coordinates  (44,35).  The  input  array  is  64x64 
pixels.  The  scene  noise  is  Gaussian,  white,  with  mean  mr,=0.5  and  standard  deviation  o=  0.3  and 
is  spatially  disjoint  with  the  target.  The  3-D  plot  of  the  output  of  the  minimum  mean  square  error 
filter  is  shown  in  Fig.4(b).  It  is  evident  that  the  minimum  mean  square  error  filter  provides  a  very 
good  result  as  evidenced  by  a  well  defined  correlation  peak  and  a  small  background  output  noise. 
This  is  because  the  desired  output  of  the  minimum  mean  square  error  filter  is  a  delta  function  which 
shaipens  the  correlation  peak  while  it  minimizes  the  output  background  noise. 
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In  summary,  we  have  described  a  minimum  mean  square  error  filter  for  pattern  recognition 
problems  with  input  scene  noise  that  is  spatially  disjoint  with  the  target.  The  filter  is  designed  to 
have  an  output  which  is  a  delta  function  located  at  the  position  of  the  target.  The  window  of  the 
reference  signal  strongly  influences  the  filter  function  and  the  filter  is  matched  to  detect  both  the 
reference  signal  and  the  spatially  disjoint  window  function  of  the  reference  signal.  The  importance 
of  the  window  information  increases  when  the  mean  of  the  scene  noise  background  is  strong. 


5  .  An  Optimal  Filter  Function  for  Spatially  Disjoint  Signal  and  Input  Scene  Noise 
Using  a  Peak  to  Noise  Ratio  Metric 

In  this  section,  an  optimum  filter  function  is  designed  for  an  input  signal  r(t)  which  is  spatially 
disjoint  from  the  input  scene  noise  h(t).  The  filter  is  designed  to  detect  a  target  r(t)  located  in  an 
unknown  random  position  x. 

We  assume  that  w,.(t)  is  the  window  function  that  defines  the  input  target.  The  window  is  unity 
where  the  function  is  defined  and  zero  otherwise,  that  is,  r(t)Wr(t)=r(t).  w„(t)  is  the  window 
function  which  defines  the  limited  size  input  scene,  and  w(t,t)  is  a  window  function  that  defines 
the  scene  noise,  that  is, 

w(t,T)=Wo(t)-Wr(t-T).  (27) 

The  Fourier  transforms  of  the  target  r(t)  and  its  window  are  denoted  by  R(co)  and  W^Cco), 
respectively.  The  Fourier  transform  of  the  input  scene  window  Wg(t)  is  denoted  by  W^(co).  The 
Fourier  transform  of  the  input  scene  noise  window  w(t,T)=Wo(t)-Wr(t-T)  is: 

W(co,T)=Wo(a))-Wj-(co)exp(-jcox)  .  (28) 
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In  deriving  the  filter  function,  the  following  assumptions  are  made: 

1)  The  scene  noise  n(t)  is  defined  as: 

n(t)=n(t)w(t,T)=n(t)[wo(t)-Wr(t-T)],  (29) 

where  n(t)  is  a  wide  sense  stationary  random  noise  with  expected  value  m^^.  The  tact  that  the  scene 
noise  is  spatially  disjoint  with  the  target  is  indicated  by  the  product  of  the  wide-sense  stationary 
noise  n(t)  and  the  window  function  w(t,T)  in  Eq.  (29).  The  scene  noise  n(t)  is  not  stationary  due  to 
the  limited  size  of  the  input  scene. 

2)  It  is  assumed  that  the  target's  random  location  "x"  has  a  uniform  probability  density  function 
f(T): 

f  1/d  xe  scene  area 

f(t)=  {  ,  (30) 

0  elsewhere 

where  d=Wo(0). 

3)  The  unknown  illumination  coefficient  a  is  a  random  variable  with  mean  a  and  variance  cl. 

4)  It  is  assumed  that  the  wide-sense  stationary  noise  n(t),  the  illumination  coefficient  a,  and  the 
target  location  x  are  statistically  independent  of  one  another. 

The  filter  impulse  response  and  transfer  function  are  denoted  by  h(t)  and  H(co),  respectively. 
The  filter  output  y(t,x)  when  the  target  is  located  at  position  t=x  in  the  input  plane  is: 

y(t,x)=|'  h(t-t')[ar(t'-x)-rn(t')w(t',x)]df  (31) 
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The  optimal  filter  function  H^pj  is  designed  to  maximize  the  ratio  of  the  expected  value  squared  of 
the  correlation  peak  to  the  average  expected  value  of  the  output  signal  energy: 


PNR  =  {E[y(x,x)]p 
E([y(t,x)]2) 


(32) 


Here,  PNR  is  the  peak-to-noise  ratio  metric,  the  symbol  in  the  denominator  "{■  |"  denotes  spatial 
averaging  (integration)  over  t,  and  E{  •}  denotes  the  statistical  expected  value. 

Given  the  target  location  x,  the  conditional  expected  value  of  the  output  signal  is: 


E[y(t,x)lx]= 


h(t-t')[  ar(t'-x)+mnW(t',x)]dt'= 


H((ja)[  aR(co)e-iw'^+mnW((jo,x)]ej"'dco, 


(33) 


where  E[.]  denotes  expected  value.  The  expected  value  of  the  correlation  peak  at  l=x  is: 


E[y(x,x)]=:E{E[y(x,x)lx]}= 


H(ca)[  aR(oo)+mnWi((o)]d(o, 


(34a) 


where 


Wi(co)=E(W(ff),x)ei“^)=|Wo(co)p/d-Wr((o),  and  d=Wo(0); 


(34b) 


The  conditional  variance  of  the  output  signal  is:!'^! 


Var[y(t,x)lx]= 


h(t-t')h(t-t")[aar(t'-x)r(t"-x)+w(t',x)w(t”,x)Cn(t'-f')]dt'dt" 
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where  Var[.]  denotes  variance,  and  Cf^(t)  and  Sj^(co)  are  the  covariance  and  the  power  spectrum  of  the 
scene  noise  n(t),  respectively. 

Let  us  denote  the  spatial  average  of  x(t)  as 

+  {E[y(t,x)lT]p  .  Provided  that  the  filter  function  is  real  and  the  spatial  extent  of  the  output  is  L,  we 
have: 


x(t)dt.  We  compute  E{[y(t,T)]^lT}  =Var[y(t,T)lT] 


where  *  denotes  convolution.  If  the  spatial  extent  of  the  output  is  equal  to  that  of  the  input,  then 
d=L. 

Using  Eqs.  (36)-(37),  we  have: 


E([y(t,T)]2|T)=V4y(t,T)lT]  +{E[y(t,T)lx]p 

[ai|R((jo)p+-iHW((i),T)|2*[Sn(co)-27tm26(to)]-fjaR((o)e-J“WmnW(co,T)pl|H((o)pdco  (38) 
2;iL  I  L  271  J 
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Given  that  t  is  uniformly  distributed,  the  expected  value  of  e( [y(t,x)]2|T}  is: 


E([y(t,T)]2)=;^  (a2|R(ro)p+;^W2(®)*Sn(Q3)+|hR(co)+mnW](0)p-m2jWi(0)pj|H((n)pd(o  (39) 


where  Wi(cd)  is  defined  in  Eq.  (34b),  and; 


W2((o)=E{|W(co,T)p)=|Wo(co)p+|Wr(co)p-2|Wo(co)pReal[Wr(®)]/d 


Therefore,  the  output  PNR  of  the  filter  is: 


PNR=- 


H((o)[5R(co)+mnW,(co)]d 


|R(®)P+^W2(co)*Sn(co)+|aR(co)+mnWi(co)p-ni2|Wi(co)p|H(®)pdco 


^  ja2|R(co)p+^ 


Applying  the  Schwarz  inequality,  we  have; 


PNR<-J- 

271 

L 


aR(co)+mnWi(co) 


a2|R((D)p+^W2(co)*Sn(co)+|aR(co)+mnWi  (co)p-m2|W|  (co)f 


The  equality  holds  when 


Hopt(®)= 


aR(co)+mnW](a)) 

a2|R(co)p+^W2(co)*Sn((D)+|aR(co)+mnWi(co)p-m2|Wi((o)p’ 


where  Hopj(oo)  is  the  complex  conjugate  of  the  filter  transfer  function. 
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Equation  (43)  describes  the  optimal  filter  function  which  detects  a  signal  in  the  presence  of 
spatially  disjoint  scene  noise.  The  filter  optimizes  the  correlation  performance  by  maximizing  the  ratio 
of  the  correlation  peak  intensity  to  the  average  output  energy.  The  filter  function  can  be  determined 
in  terms  of  the  reference  signal,  the  spatial  size  oi  the  input  scene,  the  first  and  second  order  statistics 
of  the  input  scene  noise,  and  the  input  window  functions.  The  filter  takes  into  account  the  non¬ 
whiteness  of  the  input  scene  noise,  the  non-stationarity  of  the  input  signal  due  to  the  limited  spatial 
size  of  the  scene,  the  spatial  disjointness  of  the  target  and  the  scene  noise,  and  the  possible  variation 
of  the  target  illumination. 

For  an  infinite  size  input  scene,  w(t,T)=:l-w  (t-x),  and  we  can  write: 


W,(CO)=5(CO)-Wr(CO)=Wi(CO) 


W2(co)=lWi(co)P 


In  this  case  the  optimal  filter  is: 


_ aR(m)+mnWi(co) _ 

ai|R(co)P+;:^W  I  ((o)p'f=[Sn(co)-27tm25(oo)]+|aR(co)-rmnW  i  (rajp 
2k 


In  this  case,  if  the  scene  noise  is  a  constant  [n(t)=m,J,  and  the  illumination  is  known  [a  -a  ],  then 
using  Eq.  (46),  the  optimum  filter  is: 


H(C0)=Hopt(C0)=  ^  ^ 

opiv  aR(®)+mnWi(0) 


Equation  (47)  is  an  inverse  filter  and  its  output  is  a  delta  function. 

It  can  be  seen  from  Eqs.  (43)-(46)  that  for  the  case  of  disjoint  signal  and  scene,  the  optimum  filter 
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is  veiy  different  from  the  conventional  matched  filter  expression.  The  optimum  filter  is  matched  to  the 
expected  value  of  the  input  scene  and  the  window  of  the  input  scene  noise  plays  an  important  role  in 
the  filter  function.  For  large  scene  noise,  the  scene  noise  window  function  can  be  as  important  as  the 
the  reference  function. 

If  we  ignore  the  amplitude  variation  on  the  target,  Eq.(46)  can  be  simplified  to  the  form  which  is 
exactly  the  same  as  the  MMSE  filler  in  Eq.(25).  Both  filters  attempt  to  maximize  the  output  signal  at 
the  tai'get  location  while  minimizing  the  output  signal  at  other  locations.  The  MMSE  filter  is  designed 
to  reduce  the  difference  between  ideal  della  function  at  the  target  location  and  the  filter  output.  Thus,  a 
very  well  defined  output  is  produced  at  the  target  location.  The  optimum  filter  presented  in  Eq.  (43) 
achieves  this  by  maximizing  the  output  at  the  target  location  while  minimizing  the  output  signal  energy 
at  other  locations.  Therefore,  both  approaches  result  in  similar  filter  solutions.  The  optimum  filter  is 
derived  in  continuous  form  and  takes  into  account  the  effect  of  the  limited  size  of  the  input  scene  and 
the  amplitude  variation  of  the  target.  The  MMSE  filter  is  derived  in  discrete  form. 

Computer  simulation  of  the  optimum  filter  of  Eq.(43)  is  provided  for  the  input  scene  shown  in 
Fig.  5(a).  The  tank  is  used  as  the  reference  signal.  The  target  is  placed  in  a  scene  noise  which  is 
non-stationary  due  to  the  limited  spatial  size  of  the  input[see  Fig.  5(a)].  The  input  scene  noise  n(t)  is 
white  Gaussian  noise  with  mean  of  0.5,  and  standard  deviation  of  0.3.  The  amplitude  coefficient  a  is 
assumed  to  be  one.  The  outputs  of  the  optimum  filter  [please  see  Eq.(43)],  the  conventional  matched 
filter,  and  the  phase-only  filter  for  the  same  input  scene  are  plotted  in  Figures  (5b),  (5c),  and  (5d), 
respectively. 

In  this  section,  an  optimum  filter  function  is  described  to  detect  a  target  in  the  presence  of  spatially 
disjoint  input  scene  noise.  The  filter  function  is  designed  to  optimize  the  ratio  of  the  expected  value 
squared  of  the  output  correlation  peak  to  the  average  expected  value  of  the  output  signal  energy.  This 
metric  maximizes  the  correlation  peak  and  minimizes  the  output  energy.  The  filter  function  takes 
into  account  the  scene  noise  color,  the  non-stationarity  of  the  input  scene  due  to  the  limited  spatial 
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size  of  the  input,  the  spatially  disjoint  target  and  scene  noise,  and  the  unknown  variation  of  the  target 
illumination.  The  computer  simulation  results  demonstrate  that  this  filter  function  yields  good 
performance  with  well  defined  output  correlation  peak  and  low  output  sidelobe. 


6  .  Generalization  of  the  Matched  Filter  Function  for  Spatially  Disjoint  Target  and 
Input  Scene  Noise 

Optimum  filters  such  as  matched  filters  are  designed  by  maximizing  a  metric  such  as  the 
correlation  SNR.  It  was  shown  in  section  2  that  the  classic  definition  of  the  correlation  SNR  may  not 
be  a  good  criterion  for  designing  a  filter  or  evaluating  the  performance  of  a  filter  when  the  scene  noise 
is  spatially  disjoint  with  the  target  in  the  scene.  In  this  case,  the  scene  noise  is  non-stationary  since  it 
does  not  overlap  the  target.  The  other  reason  for  the  non-stationarity  of  the  input  scene  noise  could  be 
the  limited  size  of  the  input  scene.  Thus,  a  modified  version  of  the  correlation  SNR  is  necessary  to 
design  the  filter  function. 

For  matched  filters  the  classic  definition  of  the  correlation  SNR  is  the  ratio  of  the  expected  value 
squared  of  the  correlation  peak  to  the  variance  of  the  peak.^ ' '  If  the  input  noise  is  stationary  (input 
noise  overlaps  the  target),  the  correlation  peak  variance  is  a  good  statistical  measure  of  the  variance  at 
any  other  sample  point  of  the  output  signal.  However,  statistical  measurements  of  a  non-stationary 
random  process,  such  as  mean  and  variance,  depend  on  the  sample  position.  Therefore,  when  the 
target  and  scene  noise  are  spatially  disjoint,  the  variance  of  the  correlation  peak  is  not  a  good  measure 
of  the  output  noise  fluctuations.  Thus,  to  obtain  a  measure  of  the  output  noise  fluctuations,  we  could 
consider  the  spatial  average  of  the  output  noise  variance  over  all  the  output  pixel  positions.  The 
classic  definition  of  the  correlation  SNR  is  thus  modified  to  obtain  a  new  metric  SNR,  defined  as  the 
ratio  of  the  expected  value  squared  of  the  correlation  peak  to  the  spatial  average  of  the  output  variance: 


SNR= 


|E[y(T,x)]p 

Vai[y(t,T)] 


(48) 
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The  input  signal  mode!  and  the  notations  used  in  this  section  are  the  same  as  those  established  in 
section  5.  Here,  y(t,x)  is  the  filter  output  when  the  target  is  located  at  position  T  in  the  input  plane,  the 

symbol  "(•  |"  denotes  spatial  averaging  (integration)  over  t,  and  E{  ■ }  denotes  the  statistical  expected 
value.  We  design  a  generalized  matched  filter  function  H„((jo)  which  maximizes  the  SNR  in  Eq.  (48). 
It  will  be  shown  that  the  conventional  matched  filter  designed  for  overlapping  signal  and  noise  is  a 
special  case  of  Hg(a))  when  the  input  scene  noise  is  a  wide  sense  stationary  noise  with  zero  mean  and 
is  overlapping  with  the  target. 

Using  Eq.  (36),  the  spatial  average  of  the  output  variance  is: 


Var[y(t,T)]=E(Var[y(t,i:)lx] 


f 

W  1 
'  (2Tt)2L. 


27tai|R(co)p+W2((o)*[Sn(w)-27im?i5(a))]]|H(co)pd0  (49) 


whereW2(w)  is  defined  in  Eq.(40).  Therefore,  using  Eqs.  (34a),  (48)  and  (49),  we  have; 


SNR=- 


H((x))[  aR(co)4-mnWi((j())]dq 
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L 


^  ja^|R((0)P+;^W2(G))*[Sn(C0)-2Tim2q(C0)]j|H((0)pdC0 


271 

aR(co)+mnWi(co) 


y]a2|R((0)p+^W2(C0)*[Sn((0)-27tm25((0)J[ 


d(jo 


(50) 


(51) 


where  Wi(®)  is  defined  in  Eq.(34b).  The  equality  yields  the  new  generalized  matched  filter; 


H*(co)=H*(co)= - 

a2|R(co)p+ 


aR(co)+mnWi(to) 

^W2(a))*[Sn(co)-27i:m25(co)] 


(52) 
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This  is  the  generalization  of  the  matched  filter  when  the  input  noise  does  not  overlap  the  target. 
For  an  input  scene  with  infinite  window  size,  the  generalized  matched  filter  is  [please  see  Eqs.(44) 
and  (45)]: 


Hg((0)= 


aR(co)+mnW](cD) 

a^|R(o3)P+;^W,(cn)p*[s„(co)-27rm25(m)_ 


(53) 


If  the  input  scene  noise  is  zero  mean,  and  the  illumination  is  known  (a=a=l),  we  obtain: 


H:(co)=-- - - ,  (54) 

;^Wi(CO)F*Sn(0) 

ZTI 

where  Wi(.)  is  defined  in  Eq.  (44).  Furthermore,  if  the  input  noise  is  overlapping  with  the  target, 
then  Wi((n)=5(co)  and  the  generalized  matched  filter  function  is  simplified  to: 


Hg(co)= 


R((0) 


(55) 


Equation  (55)  is  the  same  as  the  conventional  matched  filter  function  obtained  by  maximization  of 
the  classic  definition  of  SNR  under  the  condition  that  the  target  to  be  detected  is  in  the  presence  of 
zero  mean  overlapping  stationary  noise.  The  conventional  matched  filter  function  in  Eq.  (55)  is  a 
special  case  of  the  generalized  matched  filter  function  in  Eq.  (52).  When  the  stationary  input  noise  is 
overlapping  the  target,  the  variance  of  the  correlation  peak  is  a  good  statistical  measure  of  the  variance 
at  any  other  point  of  the  output  signal.  Thus,  the  spatial  average  of  the  output  variance  is  the  same  as 
the  variance  of  the  correlation  peak.  For  the  overlapping  target  and  input  noise,  the  metric  SNR  is 
equivalent  to  the  classic  definition  of  SNR,  and  the  conventional  matched  filter  solution  is  obtained 
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by  maximizing  these  metrics. 


7.  Summary 

In  many  pattern  recognition  problems,  the  input  scene  noise  is  spatially  disjoint  with  the  target  in 
the  scene.  We  have  shown  that  for  this  class  of  problems,  the  matched  filter  expressions  and  the 
optimum  receivers  derived  under  the  overlapping  input  signal  and  .scene  noise  assumption  may  not 
perform  well  in  the  presence  of  spatially  disjoint  input  signal  and  scene  noise.  The  matched  filter 
function  is  derived  by  maximizing  the  conventional  correlation  SNR  metric.  The  limitations  of  the 
conventional  SNR  metric  for  spatially  disjoint  target  and  scene  noise  are  discussed. 

In  this  report,  we  have  presented  four  possible  solutions  to  detect  targets  in  spatially  disjoint  scene 
noise.  The  first  solution  is  an  optimum  receiver  designed  based  on  multiple  alternative  hypothesis 
testing  to  detect  a  signal  in  spatially  disjoint  scene  noise.  Both  noisy  targets  (detector  noise  present) 
and  noise  free  targets  which  are  spatially  disjoint  with  the  scene  noise  and  have  an  unknown 
illumination  scale  are  discussed.  It  is  shown  that  for  a  noise  free  target,  the  actual  scene  noise 
statistics  becomes  irrelevant  to  the  detection  process.  In  this  case,  the  optimum  receiver  is  similar  to  a 
correlator  normalized  by  the  input  scene  energy  within  the  target  window. 

The  second  approach  is  based  on  Wiener  filtering.  A  minimum  mean  square  error  filter  is 
presented  for  pattern  recognition  problems  with  input  scene  noise  that  is  spatially  disjoint  with  the 
target.  The  filter  is  designed  to  have  an  output  which  is  a  delta  function  located  at  the  position  of  the 
target.  It  is  shown  that  the  filter  performs  very  well.  The  window  of  the  reference  signal  strongly 
influences  the  filter  function  and  the  filter  is  matched  to  detect  both  the  reference  signal  and  the 
spatially  disjoint  window  function  of  the  reference  signal. 

The  third  technique  is  based  on  maximizing  a  performance  metric  to  design  an  optimum  filter  to 
detect  signals  in  spatially  disjoint  scene  noise.  An  optimum  filter  function  is  designed  to  maximize  the 
ratio  of  the  correlation  peak  to  the  average  output  energy  for  a  signal  which  is  spatially  disjoint  with 


31 


the  input  scene  noise.  It  is  shown  that  the  filter  produces  a  sharp  output  signal  at  the  target  location 
with  a  low  output  noise.  Other  types  of  optimal  filters  can  be  designed  by  maximizing  variations  of 
similar  performance  metrics  which  are  defined  for  the  spatially  disjoint  target  and  scene  noise 
problem. 

The  fourth  approach  is  based  on  modifying  the  conventional  definition  of  the  correlation  peak 
SNR  to  take  into  account  the  non-overlapping  nature  of  the  target  and  the  scene  noise.  For  spatially 
disjoint  target  and  scene  noise,  the  classic  definition  of  the  correlation  SNR  is  modified  to  obtain  a 
new  metric  SNR  defined  as  the  ratio  of  the  expected  value  squared  of  the  correlation  peak  to  the 
spatial  average  of  the  output  variance.  By  maximizing  the  SNR,  we  have  designed  a  generalized 
matched  filter  function  that  can  detect  a  target  in  the  presence  of  spatially  disjoint  scene  noise.  It  is 
shown  that  the  widely  used  conventional  matched  filter  functions  designed  for  additive  signal  and 
noise  are  a  special  case  of  this  generalized  matched  filter  when  the  input  noise  is  overlapping  with  the 
target  and  has  zero  mean. 

We  wish  to  thank  Prof.  Bahaa  E.  A.  Saleh  and  V.  Laude  for  many  rewarding  discussions,  and 
A.  Fazlolahi,  G.  Zhang,  and  F.  Parchekani  for  performing  the  computer  simulation  tests. 
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Figure  Captions 

Figure  1  Input  images 

(a)  The  target  (tank)  in  the  spatially  disjoint  (non-overlapping)  additive  input  scene  noise 
used  in  the  correlation  tests. 

(b)  The  target  (tank)  in  the  overlapping  additive  input  scene  noise  used  in  the  correlation 
tests.  The  input  noise  statistics  are  the  same  as  in  Fig.  1(a). 

Figure  2  Output  correlation  intensity. 

(a)  The  correlation  output  obtained  by  a  conventional  matched  lilter  tor  the  relerence 
signal  in  the  presence  of  the  spatially  disjoint  noise  as  shown  in  Fig.  1(a). 

(b)  The  correlation  output  obtained  by  a  conventional  matched  filter  lor  the  reference 
signal  in  the  presence  of  the  additive  overlapping  noise  as  shown  in  Fig.  1  (b). 

Figure  3  Computer  simulation  results  of  an  optimum  receiver  designed  for  spatially  disjoint  target 
and  scene  noise. 

(a)  Target  (tank)  in  spatially  disjoint  (non-overlapping)  background  noise  used  in  the 
tests.  The  background  noise  is  Gaussian  with  mean  of  0.5  and  standard  deviation  of  0.3. 
The  target  has  an  illumination  of  0.5  which  is  unknown  to  the  receiver. 

(b)  Example  of  the  test  function  [-Civ[(tj)]  for  the  reference  signal  in  the  presence  of 
spatially  disjoint  background  noise  shown  in  Fig.  1(a).  C]\[(tj)  is  given  by  Eq.  (15). 

m 

(c)  Classic  correlation  receiver  solution  ^  r(ti)s(tj,-tj)  for  the  reference  signal  in  the 

i=i 

presence  of  spatially  disjoint  background  noise  shown  in  Fig.  1(a).  The  figure  illustrates 
the  pixels  around  the  correlation  peak.  The  correlation  peak  cannot  be  seen  because  it  is 
below  the  output  noise  floor. 

Figure  4  Computer  simulation  results  of  a  minimum  mean  square  error  filter  designed  for  spatially 
disjoint  target  and  scene  noise. 
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(a)  The  input  scene  with  a  target  (tank)  and  spatially  disjoint  scene  noise. 

(b)  minimum  mean  square  error  filter  output  in  response  to  the  data  in  (a). 

Figure  5  Computer  simulation  results  of  an  optimum  filter  [please  see  Eq.  (43)]  designed  for 
spatially  disjoint  target  and  scene  noise. 

(a)  Input  image  used  in  the  computer  simulations.  The  target  is  the  tank.  The  input 
scene  noise  is  white  Gaussian  with  mean  0.5,  standard  deviation  of  0.3. 

(b)  Output  of  the  optimal  matched  filter 

(c)  Correlation  output  of  the  conventional  matched  filter 

(d)  Output  of  the  phase  only  filter 
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9.  Appendix  A 

The  Effect  of  Input  Scene  Noise  Mean  on  the  Correlation  Performance  of  the 

Conventional  Matched  filter 

We  examine  the  effect  of  input  scene  noise  mean  on  the  correlation  performance  of  the 
conventional  matched  filters  for  both  the  overlapping  and  non-overlapping  (spatially  disjoint)  input 
scene.  One  dimensional  notation  is  used  for  clarity  and  simplicity. 

It  is  assumed  that  the  input  signal  contains  a  target  r(t)  and  scene  noise  n(t).  For  overlapping  input 
scene  and  signal,  we  have  n(t)=n(t),  where  n(t)  is  wide  sense  stationary.  For  non-overlapping  input 
scene  and  signal,  we  have  n(t)=[l-W|.(t)]n(t).  Here,  Wj.(t)  is  a  uniform  window  function,  that  is,  Wj.(t) 
is  unity  within  the  area  where  the  target  is  defined,  and  is  zero  elsewhere.  To  maintain  stationarity,  it 
is  assumed  that  the  input  scene  is  sufficiently  large  along  the  coordinate  axis. 

The  filter  impulse  response  is  denoted  by  h(t).  The  filter  output  is: 


h(t-t')r(t')'dt'+ 


h(t-t')n(t')dt' 


The  noisy  output  correlation  peak  is 


(a-1) 


y(0)=  h(-t')r(t')dt'+  h(-t')n(t')dt' 


(a-2) 


The  difference  between  the  correlation  peak  and  the  output  signal  at  any  position  "t"  is: 


[h(-t')-h(t-t')]r(t')dt'+  [h(-t')-h(t-t')]H(t')dt' 


(a-3) 


37 


We  consider  the  difference  As  to  investigate  the  filter  performance  for  detecting  the  target  r(t)  in  the 
presence  of  the  scene  noise  n(t).  Obviously,  a  large  As  indicates  that  the  peak  is  much  higher  than  the 
output  sidelobe. 

Let  us  assume  that  the  input  noise  n(t)  is  white  with  mean  mj^.  We  examine  the  effect  of  the  noise 
mean  on  the  correlation  performance  of  the  conventional  matched  filter,  that  is  h(t)=r(-t),  for  both  the 
the  overlapping  and  the  non-overlapping  input  signal  and  scene  noise  cases.  The  expected  value  of 
the  output  correlation  peak  using  Eq.  (a-2)  is: 


^yO''(0)]=R?r(0)+mn  r(f)dt' 

(a-4) 

J 

E[ynov(0)]=RC.(0) 

ta-5) 

The  superscripts  "ov"  and  "nov"  denote  the  output  measurements  corresponding  to  the  overlapping 
and  non-overlapping  input  cases,  respectively.  The  expected  value  of  the  difference  As  between  the 
output  correlation  peak  and  the  output  correlation  is: 


E[As°1=R?r(0)-R?r(t)+mnj  [r(t')-r(-t+t')]dt'-R^r(0)-R?r(t) 

E[As"°1=R|?r(0)-R|;r(t)-mn  f  r(-t+f)[l-Wr(t)]dt' 


(a-6) 

(a-7) 


The  second  equality  in  Eq.  (a-6)  holds  under  the  assumption  that  the  bounds  of  the  integral  are 
sufficiently  large. 

For  large  values  of  t  in  Eqs.  (a-6)-(a-7),  r(-t+t')  is  shifted  out  of  the  window  W|,(f),  and  we  have: 


E[As°1=Rf,(0) 


(a-8) 


E[As"°1=R?r(0) 


r(f  )df  <  e[as°'' 


(a-9) 


38 


Here,  it  is  assumed  that  R,.,.(t)  approaches  zero  for  large  t. 

It  can  be  seen  from  Eqs.  (a-4),  (a-6),  and  (a-8)  that  for  the  overlapping  signal  and  scene  noise 
case,  the  output  correlation  peak  is  affected  by  the  noise  mean,  and  the  difference  between  the  output 
correlation  peak  and  the  noisy  output  sidelobe,  As,  is  independent  of  the  noise  mean. 

In  contrast,  it  can  be  seen  from  Eq.  (a-5)  that  for  the  non-overlapping  (spatially  disjoint)  signal 
and  scene  noise  case,  the  correlation  peak  of  the  conventional  matched  filter  is  independent  of  the 
noise  mean.  However,  it  can  be  seen  from  Eqs.  (a-7)  and  (a-9)  that  the  difference  between  the 
correlation  peak  and  the  noisy  output  sidelobe.  As,  is  affected  by  the  noise  mean.  Inequality  (a-9) 
shows  that  the  correlation  performance  of  the  conventional  matched  filter  for  the  spatially  disjoint 
input  signal  and  scene  noise  is  inferior  to  that  of  the  overlapping  input  case.  In  addition,  the  filter 
performance  for  the  spatially  disjoint  input  signal  and  scene  noise  deteriorates  as  the  input  scene  noise 
mean  value  increases. 

The  correlation  performance  of  the  conventional  matched  filter  for  the  spatially  disjoint  input 
signal  and  scene  noise  case  is  improved  when  the  mean  of  the  scene  noise  is  small.  For  instance, 
when  the  scene  noise  mean  is  zero,  e[As""'']  is  maximized.  However,  real  images  have  a  non-zero 
mean.  In  a  filter  based  correlation  system,  using  a  DC  block  in  the  Fourier  plane  can  remove  the 
mean  of  the  input  signal.  Under  this  condition,  each  output  signal  pixel  value  is  shifted  down  by  the 
same  amount.  This  amount  is  determined  by  the  mean  of  the  input  signal  and  the  mean  of  the 
reference  signal.  However,  the  expected  value  of  the  difference  between  the  correlation  peak  and  the 
output.  As,  is  not  affected.  Furthermore,  due  to  the  spatial  size  of  DC  block,  this  method  will 
remove  the  low  frequency  information  of  the  input  scene.  If  the  target  has  significant  energy  in  the 
low  spatial  frequencies,  a  DC  block  can  result  in  loss  of  light  efficiency,  and  loss  of  information 
about  the  target. 
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10.  Appendix  B 


List  of  Symbols 


a 

a 

/N 

a 

Ci(co) 

Cn(t) 

E{-} 

f(t) 

h(t) 

Ha,S) 

H(co) 

Hg(co) 

Hki(®) 


Hopt(co) 

k 

m 

MSE 

n(ti) 

n(t) 


unknown  illumination  coefficient  of  the  target  in  the  input  scene 
mean  of  the  random  variable  a 
maximum  likelihood  estimate  of  a 

Fourier  component  which  generates  the  first  order  correlation  output 
covariance  of  the  input  noise  n(t) 
expected  value 

probability  density  function  of  the  target  location  in  the  input  scene 
impulse  response  of  the  filter 

hypothesis  H(j,a)  correspond  to  the  target  located  at  position  tj  and  with 
amplitude  a 

Fourier  transform  of  h(t),  filter  transfer  function 
transfer  function  of  the  generalized  matched  filter 

first  order  component  of  the  kth  law  nonlinearly  transformed  matched  filter 
function 

optimal  filter  transfer  function 

degree  of  the  kth  law  nonlinearity 

number  of  the  samples  taken  over  the  input  scene  s(t) 

expected  value  of  the  input  noise  n(t) 

number  of  the  samples  taken  over  the  target  r(t) 

mean  square  error 

sample  of  the  noise  n(t) 

wide  sense  stationary  noise  in  the  input  plane 
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nB(0 

nd(t) 

iijCt) 

no(t) 

H(t) 

N((£)) 

N(03) 

Pnj(') 

Pnd^'^ 

P(') 

PNR 

^•k(0 

r(t) 

r(ti) 

R(co) 

R(v) 

R?r(ti) 

R^r(ti)  or  Rf^, 
R"nn(t,) 

S(tj) 

S(t) 


spatial  disjoint  background  noise  in  the  input  plane 
detector  noise 

spatially  disjoint  input  scene  noise  in  section  3  which  is  the  sum  of  the 

background  noise  and  the  detector  noise 

input  noise  minus  by  its  mean  value 

disjoint  input  scene  noise  in  section  5,6  and  appendix  A 

output  noise  contribution  to  the  output  correlation  peak 

Fourier  transform  of  the  noise  n(t) 

Fourier  transform  of  the  spatially  disjoint  input  scene  noise  n(t) 
joint  probability  density  function  of  samples  of  nj(t) 
joint  probability  density  function  of  samples  of  11^(1) 
likelihood  function 

the  peak  to  noise  ratio  defined  as  the  ratio  of  the  expected  value  squared  of  the 

correlation  peak  to  the  average  expected  value  of  the  output  signal  energy 

inverse  Fourier  transform  of  |R(co)pexp{j(t)R(to)} 

target  function  in  the  input  scene 

sample  of  the  target  function 

Fourier  transform  of  the  target  function  r(t) 

discrete  Fourier  transform  of  the  target  sample  sequence  [r(tj)] 

correlation  of  the  sample  sequence  of  the  target  [r(tj)] 

(ti)  correlation  between  [r(tj)]  and  [w(tj)] 
correlation  of  the  sample  sequence  of  the  noise  [n(tj)] 
sample  of  the  input  scene 
input  scene  function 

m^.  samples  taken  over  the  target  in  the  input  scene 
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sj  (m-m^.)  samples  taken  over  the  input  scene  noise 

S((x))  Fourier  transform  of  the  input  scene  function  s(t) 

S(v)  discrete  Fourier  transform  of  the  input  scene  sample  sequence  [s(t-)] 

Sj^(co)  power  spectrum  of  the  noise  n(t) 

SNR  conventional  definition  of  the  signal  to  noise  ratio  defined  as  the  ratio  of  the 

expected  value  squared  of  the  correlation  peak  to  the  variance  of  the  peak 
SNR  modified  signal  to  noise  ratio  for  spatially  disjoint  target  and  scene  noise 

defined  as  the  ratio  of  the  expected  value  squared  of  the  correlation  peak  to  the 
spatial  average  of  the  output  variance 

W|.(t)  window  function  which  defines  the  target  in  the  input  scene 

WQ(t)  window  function  which  defines  the  input  scene 

w(t)  window  function  which  defines  the  spatially  disjoint  noise 

w(tj)  sample  of  the  window  function  which  defines  the  spatially  disjoint  input  noise 

W((o)  Fourier  transform  of  w(t) 

W[.(co)  Fourier  transform  of  W|.(t) 

Wo(ro)  Fourier  transform  of  Wo(t) 

y(t;)  discrete  form  of  ideal  filter  output  in  response  to  an  input  with  spatially  disjoint 

signal  and  scene  noise 

y(ti)  discrete  form  of  filter  output  in  response  to  an  input  with  spatially  disjoint  signal 

and  scene  noise 

y(t)  filter  output  in  response  to  an  input  with  spatially  disjoint  target  and  scene  noise 

As  difference  between  the  correlation  peak  and  the  output  signal  at  any  position 

(t)N(©)  Fourier  phase  of  the  spatially  disjoint  input  noise  n(t) 

(1)r(co)  Fourier  phase  of  the  target  r(t) 

(})s(co)  Fourier  phase  of  the  input  scene  s(t) 

On  standard  deviation  of  the  noise  n(t) 
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(a) 


(b) 

Figure  1  Input  images. 

(a)  The  target  (tank)  in  the  spatially  disjoint  (nonoverlapping)  additive  input  scene 
noise  used  in  the  correlation  tests. 

(b)  The  target  (tank)  in  the  overlapping  additive  input  scene  noise  used  in  the 
correlation  tests.  The  input  noise  statistics  are  the  same  as  in  Fig.  1(a). 
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Figure  2  Output  correlation  intensity 

(a)  The  correlation  output  obtained  by  a  conventional  matched  filter  for  the  reference 
signal  in  the  presence  of  the  additive  spatially  disjoint  noise  as  shown  in  Fig.  1(a) 

(b)  The  correlation  output  obtained  by  a  conventional  matched  filter  for  the 
reference  signal  in  tlte  presence  of  the  additive  overlapping  noise  as  shown  in 


Figure  3  Computer  simulation  results  of  an  optimum  receiver  designed  for  spatially  disjoint  target 


and  scene  noise, 

(a)  Target  (tank)  in  spatially  disjoint  (non-overlapping)  background  noise  used  in  the 
tests.  The  background  noise  is  Gaussian  with  mean  of  0.5  and  standard  deviation  of  0.3. 
The  target  has  an  illumination  of  0.5  which  is  unknown  to  the  receiver. 

(b)  Example  of  the  test  function  [-Ci\i(tj)]  for  the  reference  signal  in  the  presence  of 
spatially  disjoint  background  noise  shown  in  Fig.  1(a).  C[y(tj)  is  given  by  Eq.  (15). 

ni 

(c)  Classic  correlation  receiver  solution  ^  r(tj)s(ti-tj)  for  the  reference  signal  in  the 

i=i 

presence  of  spatially  disjoint  background  noise  shown  in  Fig.  1(a).  The  figure  illustrates 
the  pixels  around  the  correlation  peak.  The  correlation  peak  cannot  be  seen  because  it  is 
below  the  output  noise  floor. 
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(b) 

Figure  4  Computer  simulation  results  of  a  minimum  mean  square  error  filter  designed  for  spatially 
disjoint  target  and  scene  noise. 

(a)  The  input  scene  with  a  target  (tank)  and  spatially  disjoint  .scene  noise. 

(b)  minimum  mean  square  error  filter  output  in  response  to  the  data  in  (a). 


47 


A 

0.1 

0  , 

.  5 

O.l 

1  - 

•J 

0 

0 

Figure  5  Computer  simulation  results  of  an  optimum  filter  [please  see  Eq.  (43)]  designed  for 
spatially  disjoint  target  and  scene  noise. 

(a)  Input  image  used  in  the  computer  simulations.  The  target  is  the  tank.  The  input 
scene  noise  is  white  Gaussian  with  mean  0.5,  standard  deviation  of  0.3. 

(b)  Output  of  the  optimal  matched  filter 

(c)  Correlation  output  of  the  conventional  matched  filter 

(d)  Output  of  the  phase  only  filter 
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